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Abstract. From a detailed elemental abundance analysis of 102 F and G dwarf stars we present abundance trends in the 
Galactic thin and thick disks for 14 elements (O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Zn, Y, Ba, and Eu). Stellar parameters 
and elemental abundances (except for Y, Ba and Eu) for 66 of the 102 stars were presented in our previous studies (Bensby et 
al. 2003"2004a^. The 36 stars that are new in this study extend and confirm our previous results and allow us to draw further 
conclusions regarding abundance trends. The s-process elements Y and Ba, and the r-element Eu have also been considered 
here for the whole sample for the first time. With this new larger sample we now have the following results: 1) Smooth and 
distinct abundance trends that for the thin and thick disks are clearly separated; 2) The Q--element trends for the thick disk show 
typical signatures from the enrichment of SN la; 3) The thick disk stellar sample is in the mean older than the thin disk stellar 
sample; 4) The thick disk abundance trends are invariant with galactocentric radii (/?,„); 5) The thick disk abundance trends 
appear to be invariant with vertical distance (Z^jnx) from the Galactic plane. Adding further evidence from the literaure we argue 
that a merger/interacting scenario with a companion galaxy to produce a kinematical heating of the stars (that make up today's 
thick disk) in a pre-existing old thin disk is the most likely formation scenario for the Galactic thick disk. 
The 102 stars have -1 < [Fe/H] < -1-0.4 and are all in the solar neighbourhood. Based on their kinematics they have been 
divided into a thin disk sample and a thick disk sample consisting of 60 and 38 stars, respectively. The remaining 4 stars have 
kinematics that make them kinematically intermediate to the two disks. Their chemical abundances also place them in between 
the two disks. Which of the two disk populations these 4 stars belong to, or if they form a distinct population of their own, can 
at the moment not be settled. The 66 stars from our previous studies were observed with the FEROS spectrograph on the ESO 
1.5-m telescope and the CES spectrograph on the ESO 3.6-m telescope. Of the 36 new stars presented here 30 were observed 
with the SOFIN spectrograph on the Nordic Optical Telescope on La Palma, 3 with the UVES spectrograph on VLT/UT2, and 
3 with the FEROS spectrograph on the ESO 1.5-m telescope. All spectra have high signal-to-noise ratios (typically S/N > 250) 
and high resolution (R ~ 80 000, 45 000, and 1 10 000 for the SOFIN, FEROS, and UVES spectra, respectively). 

Key words. Stars: fundamental parameters - Stars: abundances - Galaxy: disk - Galaxy: formation - Galaxy: abundances - 
Galaxy: kinematics and dynamics 



1. Introduction 

During the last few years several studies have used detailed 
abundance analysis in order to establish the chemical proper- 
ties of the thick disk stellar population (e.g., Bensby et al. i2003. 
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* Based on observations collected at the Nordic Optical Telescope 
on La Palma, Spain, and at the European Southern Observatories on 
La Silla and Paranal, Chile, Proposals #65.L-0019(B), 67.B-0108(B), 
69.B-0277. 

** Tables|31Hl and|5|are only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat?J/A+A/ 



I2004al Feltzing et al. l2003l Reddy et al l2003l Tautvaisiene et 
al. I2OOTI Mashonkina & Gehren "200?; Gratton et al. "20001 
Prochaska et al. 2000; Chen et al. 2000 Fuhrmann 1998 1. 
Although the various studies take different approaches to defin- 
ing the stellar samples and though some of them are only con- 
cerned with one of the disks, there is a general agreement on 
the following: 1) The thick disk is, at a given [Fe/H], more en- 
hanced in the a-elements than the thin disk; 2) The abundance 
trend in the thin disk is a gentle slope, and 3) The solar neigh- 
bourhood thick disk stars that have been studied so far are all 
old. 

The aim of the present study is to verify and extend these re- 
sults and to add new elements into the discussion; the r-process 
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Fig. 1. Toomre diagram for the full stellar sample (102 stars). 
Thick and thin disk stars are marked by filled and open sym- 
bols, respectively. Stars that have been observed with SOFIN 
or UVES are marked by triangles and those from Bensby 
et al. <2003l l are marked by circles. "Transition objects" are 
marked by "open stars". 

element europium (Eu) and the two i-process elements yttrium 
(Y) and barium (Ba). By studying Eu and Ba Mashonkina & 
Gehren (2001 1 found that AGB stars have contributed to the 
chemical enrichment of the thick disk. By including these ele- 
ments we will be able to confirm this important finding. These 
elements will also be combined with the a-elements, in par- 
ticular our oxygen abundances from Bensby et al. ( 2004aj) . to 
shed new light on the chemical enrichment. 

The paper is organized as follows. In Sect. |2lwe describe 
the stellar sample. Section|3]describes the observations and the 
data reductions. Sections |4] and [S] briefly describe the stellar 
model atmospheres and the elemental abundance determina- 
tion. The interested reader is refered to Bensby et al. i200 3 1 for 
a detailed discussion. Section|S]describes how we determined 
stellar ages. The resulting elemental abundance trends are then 
presented in Sect.0and combined with the results from Bensby 
et al. (120031 12004at for an extended discussion. Conclusions 
and a final summary are given in in Sect. |8] The paper ends 
with an Appendix that includes a discussion of the assumptions 
made about the parameters that are used in the kinematical se- 
lection criteria for the stellar samples. 

2. Stellar sample 

Given that we, in principle, never can select a sample of local 
thick disk stars that is guaranteed to be completely free from 
intervening thin disk stars, we argue that we should keep the se- 
lection criteria as simple and as transparent as possible. In this 
sense the simplest and most honest selection is based only on 



Table 1. Kinematical data for the new stellar sample. Column 1 
gives the Hipparcos number; col. 2 gives the radial velocities 
(as measured by us); cols. 3-5 give the space velocities relative 
to the local standard of rest; cols. 6 and 7 give the calculated 
TD/D and TD/H ratios. The radial velocity for HIP 116740 
have been taken from Barbier-Brossat et al. (I1994>. 



HIP 



Ulsr Vlsr 
— [kms-'] - 



Wlsr TD/D TD/H 







THIN DISK STARS 






699 


-14.9 


37.6 


-4.9 


-0.8 


0.013 


>999 


910 


15.0 


29.0 


-7.8 


-11.1 


0.014 


>999 


2235 


-9.9 


23.2 


-26.2 


-6.7 


0.018 


>999 


2787 


11.7 


-16.5 


-21.7 


-12.9 


0.017 


>999 


3909 


7.6 


31.4 


2.5 


-4.8 


0.013 


>999 


5862 


14.3 


-34.1 


-17.1 


-6.0 


0.016 


>999 


10306 


5.9 


-10.3 


-6.8 


11.2 


0.011 


>999 


15131 


6.7 


3.6 


20.2 


-9.6 


0.017 


>999 


18833 


-29.5 


50.8 


3.6 


-14.0 


0.027 


>999 


88945 


-14.4 


-4.3 


1.1 


-3.3 


0.009 


>999 


92270 


3.6 


12.0 


9.0 


4.4 


0.011 


>999 


93185 


-40.4 


-29.0 


-19.1 


2.2 


0.014 


>999 


96258 


1.0 


29.3 


7.3 


-5.4 


0.013 


>999 


107975 


19.3 


23.6 


21.3 


-0.3 


0.018 


>999 


113174 


-28.1 


-1.8 


-25.4 


10.9 


0.017 


>999 






THICK DISK STARS 






11309 


-8.5 


73.3 


-34.2 


-55.3 


11.7 


>999 


12306 


-99.7 


146.0 


-25.7 


-33.9 


46.1 


819 


15510 


90.3 


-69.4 


-89.3 


-24.0 


73.7 


765 


16788 


-29.1 


78.0 


-2.2 


-50.3 


2.7 


>999 


18235 


100.7 


-17.0 


-158.1 


-20.1 


>999 


13.4 


20242 


-31.6 


90.5 


-59.7 


-19.8 


3.3 


>999 


21832 


111.5 


-124.6 


-44.6 


18.9 


7.3 


>999 


26828 


75.8 


-61.4 


-81.2 


40.2 


80.9 


838 


36874 


-134.9 


149.7 


-31.1 


-40.7 


179 


584 


37789 


-5.4 


4.5 


-103.4 


-15.8 


115 


657 


40613 


113.4 


-28.4 


-139.7 


-35.9 


>999 


41.0 


44075 


119.5 


-38.2 


-85.2 


76.8 


>999 


198 


44860 


66.8 


-58.9 


-68.7 


-56.7 


173 


796 


112151 


-8.3 


-10.9 


-65.9 


-46.9 


9.6 


>999 


116421 


-112.1 


-55.1 


-107.9 


69.8 


>999 


78.2 


116740 


-32.0 


-68.5 


-27.8 


47.3 


1.94 


>999 


118010 


3.3 


-45.9 


-15.8 


62.5 


8.0 


>999 




"TRANSITION OBJECTS" 






3170 


1.9 


-102.8 


-12.2 


-27.6 


0.71 


>999 


44441 


72.8 


-110.5 


-4.2 


-20.2 


0.62 


>999 


95447 


-101.3 


-107.1 


-25.7 


-13.3 


0.55 


>999 


100412 


17.0 


46.1 


-43.0 


41.8 


1.04 


>999 



the kinematics of the stars. This is also the least model depen- 
dent method. The selection method we used is described and 
discussed in Bensby et al. (.2003J . see also Bensby et al. (I2004al 
I2()()4b> . 

Our study contains two major stellar samples. They have 
been defined to kinematically resemble the thin and thick 
Galactic disks, respectively. As mentioned the criteria and 
method are described in Bensby et al. (2003). However, while 
we then used a 6 % for the normalization of the thick disk stars 
in the solar neighbourhood we here use 10 % (and consequently 
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Table 2. Wavelength coverage for the different spectral orders 
(SO) for the two settings (Blue and Red) of the CCD. 



so 


Blue 


Red 


so 


Blue 


Red 




[A] 


[A] 




[A] 


[A] 


50 






4530 


4567 


37 


6030 


-6084 


6121 - 


6172 


49 


4552 


4594 


4622 


4660 


36 


6198 


- 6253 


6291 - 


6344 


48 


4648 


4690 


4718 


4758 


35 


6375 


-6430 


6471 - 


6524 


47 


4748 


4789 


4819 


4858 


34 


6562 


-6620 


6660- 


6716 


46 


4850 


4893 


4924 


4964 


33 


6760 


-6820 


6863 - 


6920 


45 


4958 


5002 


5035 


5074 


32 


6972 


-7035 


7077- 


7136 


44 


5072 


5116 


5148 


5190 


31 


7197 


-7260 


7305 - 


7366 


43 


5190 


5235 


5267 


5310 


30 


7437 


-7503 


7550- 


7612 


42 


5312 


5360 


5392 


5437 


29 


7695 


-7762 


7810- 


7875 


41 


5442 


5490 


5526 


5570 


28 


7968 


-8040 


8088- 


8155 


40 


5578 


5628 


5662 


5709 


27 


8265 


- 8337 


8388 - 


8457 


39 


5722 


5772 


5807 


5855 


26 


8582 


- 8657 


8710- 


8783 


38 


5872 


5924 


5960 


6010 













the normalization for the thin disk population is 90 %). Reasons 
for this are given in the Appendix. 

Our total stellar sample contains 38 thick disk stars and 60 
thin disk stars. The new sample contains 17 thick disk stars, 15 
thin disk stars (all with [Fe/H] < 0) and a further 4 stars with 
kinematics intermediate between the thin and thick disks. By 
intermediate we mean that they can be classified either as thin 
disk or as thick disk stars depending on the value for the solar 
neighbourhood thick disk stellar density. Using a value of 6 % 
will classify them as thin disk stars whereas a value of 14 % 
will classify them as thick disk stars. Due to this ambiguity we 
will treat these stars (HIP 3170, HIP 44441, HIP 95447, and 
HIP 100412) separately from the two other samples and label 
them as "transition objects". The remaining 21 thick disk and 
45 thin disk stars were analyzed in Bensby et al. (I2003t . 

Radial velocities were determined for 35 of the 36 stars 
in the new sample. Good agreement to the radial velocities in 
the compilation by Barbier-Brossat et al. ( 1994i (which were 
used when selecting the stars for the observations) is generally 
found, with the exception of HIP 18833 where our radial ve- 
locity is 13kms ' lower. The average difference for the other 
stars is +0.4±3kms"^ with our measurements giving the larger 
values. For one star (HIP 1 16740) we adopted the radial veloc- 
ity as given in Barbier-Brossat et al. ( 1994 1 since there was an 
offset in the wavelength shift between the red and the blue set- 
tings in the SOFIN spectra for this star (compare Table|2j- All 
kinematical properties and the calculated TD/D and TD/H ra- 
tios (using the 10 % normalization, see Appendix) for the new 
sample, are given in Tabled 

3. Observations and data reductions 

3.1. SOFIN data 

Observations were carried out with the Nordic Optical 
Telescope (NOT) on La Palma, Spain, during two observ- 
ing runs in August (3 nights) and November (5 nights) 2002. 
The SOFIN (SOviet FINnish) spectrograph was used to ob- 
tain spectra with high resolving power (R ~ 80 000) and high 
signal-to-noise ratios (S /N > 250). A solar spectrum was also 
obtained by observing the Moon. To avoid long exposure times 



and thus the effects of cosmic rays the exposures were split into 
two or three exposures (not longer than ~ 20 min). 

The spectra were reduced using the 4A package 
(Ilyin 2000 1. This comprises a standard procedure for data re- 
duction and includes bias subtraction, estimation of the vari- 
ances of the pixel intensities, correction for the master flat 
field, scattered light subtraction with the aid of 2D-smoothing 
splines, definition of the spectral orders, and weighted inte- 
gration of the intensity with elimination of cosmic rays. The 
wavelength calibration was done using ThAr comparison spec- 
tra, one taken before and one after each individual object ex- 
posure. A typical eiTor of the ThAr wavelength calibration is 
about 10ms"' in the image center. 

In order to get large enough spectral wavelength coverage 
we observed each star twice with different settings for the CCD 
(see Table|2j. Each setting resulted in ~ 45 spectral orders. In 
this study we use spectral lines with wavelengths ranging from 
~ 4500 A to ~ 8800 A, i.e., spectral orders 26-50. For each 
star we analyzed approximately 260 spectral lines, which form 
a sub-set of the 450 lines that were analyzed for each star in 
Bensby et al. (l2003t . 

In total we observed 41 stars wifli NOT/SOFIN, but unfor- 
tunately we had to reject 1 1 of them from the analysis because 
their rotational velocities (v ■ sin /) were too high to allow equiv- 
alent width measurements (HIP 3641, HIP 4989, HIP 5034, 
HIP 6669, HIP 6706, HIP 18859, HIP 24109, HIP 45879, and 
HIP 87958) or because they were found to be spectroscopic 
binaries (HIP 17732 and HIP 109652). 

3.2. FEROS and UVES data 

Spectra for 69 stars (abundances for 66 of these were presented 
in Bensby et al. i2003t the remaining three will be discussed 
here and are labeled as "transition objects") were obtained 
with the FEROS spectrograph on the ESO 1.52-m telescope 
on La Silla in Chile in September 2000 and August/September 

2001. These spectra have R ~ 48 000 and S /N ~ 150-250 
with a wavelength coverage that is complete from ~ 4000 A to 
~ 9400 A. In each stellar spectrum we analyze a total of ~ 450 
spectral lines. The reductions and the analysis of these stars 
were presented in Bensby et al. <2003l l. 

Spectra for three additional stars were obtained with the 
UVES spectrograph on the VLT/Keuyen 8-m telescope in July 

2002. The spectra have R ~ 110000 and S /N > 350. The re- 
ductions of these observations will be discussed in a forthcom- 
ing paper where the majority of the stars (bulge and thick disk 
in situ giants) from that observing run will be presented. The 
setting of the CCD gives a wavelength coverage from ~ 5540 A 
to ~ 7560 A (with a gap between 6520-6670 A). This resulted 
in that ~ 200 spectral lines were analyzed (again a subset of the 
450 lines analyzed in Bensby et al. 2003^ . 

4. Stellar model atmospheres 

The calculation of stellar atmospheric models were done with 
the Uppsala MARCS code (Gustafsson et al. 1975 ; Edvardsson 
et al. .l993 ; Asplund et al. .1997,) . The iterative process to tune 
the stellar parameters used in the construction of the model 
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Table 3. Our program stars. Columns 1-3 give the identifications for each star, Hipparcos, HD, and HR numbers; col. 4 gives 
the spectral class as listed in the SIMBAD database; cols. 5-7 give V magnitude, parallax (n), and accuracy of the parallax {o-„), 
all from the Hipparcos catalogue; cols. 8-10 give the stellar atmospheric parameters, metallicity ([Fe/H]), effective temperature 
(Teff), and surface gravity (logg); col. 1 1 gives the microturbulence (^t); col. 12 gives the stellar mass (Al); col. 13 the bolometric 
correction (BV). The last column indicates which instrument was used to obtain the spectrum. 



Identifications 




Sp. type 


V 


n 




[Fe/H] 








M 


BC 


Spec. 


HIP 


HD 


HR 






[mas] 


[mas] 




[Kl 




[kms"'] 


[Mnl 
















THIN DISK STARS 














699 


400 


17 


F8IV 


6.21 


30.26 


0.69 


-0.20 


6250 


4.19 


1.35 


1.27 


-0.09 


SOFIN 


910 


693 




r J V 




52.94 


0.77 


-0.36 


ozzu 




1 A'^ 


1 1 n 






2235 


2454 


107 


F6V 


6.05 


27.51 


0.86 


-0.28 


6645 


4.17 


1.75 


1.30 


-0.07 


SOFIN 


2787 


3229 


143 


F5IV 


5.94 


17.97 


0.74 


-0.11 


6620 


3.86 


1.70 


1.68 


-0.06 


SOFIN 


3909 


4813 


235 


F7IV 


5.17 


64.69 


1.03 


-0.06 


6270 


4.41 


1.12 


1.17 


-0.08 


SOFIN 


5862 


7570 


370 


F8V 


4.97 


66.43 


0.64 


0.17 


6100 


4.26 


1.10 


1.04 


-0.08 


UVES 


10306 


13555 


646 


F5V 


5.23 


33.19 


0.85 


-0.17 


6560 


4.04 


1.75 


1.48 


-0.07 


SOFIN 


15131 


20407 




GIV 


6.75 


41.05 


0.59 


-0.52 


5834 


4.35 


1.00 


0.85 


-0.15 


UVES 


18833 


25322 




F5V 


7.82 


11.47 


1.02 


-0.52 


6370 


3.99 


1.75 


1.19 


-0.10 


SOFIN 


88945 


166435 




GO 


6.84 


39.62 


0.68 


-0.05 


5690 


4.37 


1.33 


0.95 


-0.15 


SOFIN 


92270 


174160 


7079 


F8V 


6.19 


34.85 


0.68 


-0.06 


6370 


4.32 


1.50 


1.20 


-0.07 


SOFIN 


93185 


176377 




GO 


6.80 


42.68 


0.64 


-0.28 


5810 


4.40 


0.90 


0.84 


-0.14 


SOFIN 


96258 


184960 


7451 


F7V 


5.71 


39.08 


0.47 


-0.02 


6380 


4.25 


1.48 


1.25 


-0.08 


SOFIN 


107975 


207978 


8354 


F6IV 


5.52 


36.15 


0.69 


-0.53 


6460 


4.06 


1.50 


1.10 


-0.10 


SOFIN 


113174 


216756 


8718 


F5II 


5.91 


24.24 


0.68 


-0.11 


6870 


4.14 


1.95 


1.51 


-0.05 


SOFIN 












THICK DISK STARS 














1 1 ^no 

1 i OUy 






F5 


7.36 


15.05 


0.91 


-0.32 


6210 


3.98 


1.40 


1.12 


-0.10 


SOFIN 


12306 


16397 




GOV 


7.36 


27.89 


1.12 


-0.53 


5765 


4.20 


0.90 


0.78 


-0.16 


SOFIN 


15510 


20794 


1008 


G8V 


4.26 


165.02 


0.55 


-0.41 


5480 


4.43 


0.75 


0.82 


-0.20 


UVES 


16788 


22309 




nn 
ou 


7 f\'\ 
1 .DJ 


22.25 


1.14 


-0.32 


5920 




1 no 


yj.oy 


— U. LJ 




18235 


24616 




G8IV 


6.68 


15.87 


0.81 


-0.71 


5000 


3.13 


0.95 


0.79 


-0.33 


SOFIN 


20242 


27485 




GO 




14.79 


0.98 


-0.26 


5650 


3.94 


1.00 


1.02 


-0.16 


SOFIN 


21832 


29587 




G2V 


7.29 


35.31 


1.07 


-0.61 


5570 


4.27 


0.65 


0.72 


-0.19 


SOFIN 


26828 


37739 




F5 


7.92 


12.32 


0.99 


-0.37 


6410 


4.15 


1.50 


1.31 


-0.09 


SOFIN 


36874 


60298 




G2V 


7.37 


25.42 


0.99 


-0.07 


5730 


4.22 


0.90 


0.98 


-0.14 


SOFIN 


37789 


62301 




F8V 


6.74 


29.22 


0.96 


-0.67 


5900 


4.09 


1.20 


0.88 


-0.15 


SOFIN 


40613 


69611 




F8 


7.74 


20.46 


1.16 


-0.63 


5740 


4.11 


0.92 


0.84 


-0.16 


SOFIN 


44075 


76932 


3578 


F8IV 


5.80 


46.90 


0.97 


-0.91 


5875 


4.10 


1.10 


0.86 


-0.17 


SOFIN 


44860 


78558 




G3V 


7.29 


27.27 


0.91 


-0.45 


5690 


4.19 


0.82 


0.90 


-0.16 


SOFIN 


112151 


215110 




G5 


7.73 


11.30 


0.93 


-0.42 


5035 


3.43 


0.85 


1.15 


-0.32 


SOFIN 


116421 


221830 




F9V 


6.86 


30.93 


0.73 


-0.45 


5700 


4.15 


0.95 


0.93 


-0.16 


SOFIN 


116740 


222317 




G2V 


7.04 


20.27 


0.76 


0.05 


5740 


3.96 


1.15 


1.12 


-0.13 


SOFIN 


118010 


224233 




GO 


7.67 


20.01 


0.74 


-0.07 


5795 


4.17 


1.00 


1.03 


-0.13 


SOFIN 












"TRANSITION OBJECTS" 












3170 


3823 


176 


GIV 


5.89 


39.26 


0.56 


-0.34 


5970 


4.11 


1.40 


1.05 


-0.12 


FEROS 


44441 


77408 




F6IV 


7.03 


19.85 


0.86 


-0.28 
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atmospheres and the abundance analysis is fully described in 
Bensby et al. J2003> . In summary the main ingredients are: 
Effective temperature (T^ff) is determined by requiring Fez 
lines with different lower excitation potentials to give equal 
abundances; Stellar mass (M) is estimated from the evolu- 
tionary tracks by Yi et al. (2003i. These two parameters are 
then used together with parallaxes and magnitudes from the 
Hipparcos catalogue (ESA I1997t to determine the surface grav- 
ity (logg) of the star (see Eq. 4 in Bensby et al. l2003> . The bolo- 



metric correction (BC) which also is needed to determine log g 
is found by interpolating in the grids by Alonso et al. J1995> . 
The microturbulence parameter (^t) is determined by forcing all 
Fe I lines to give the same abundance regardless of line strength, 
i.e. log(W,/^). 

Since the Fe ii lines have not been used in the determina- 
tion of the atmospheric parameters they can be used to check 
the derived parameters as well as the Fei abundances. This is 
an important test since the Fe i abundancs can be affected by 
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Fig. 2. [Fei/Feii] versus [Fe/H], log Tetf , \ogg, and Thin 
and thick disk stars are marked by open and filled circles, re- 
spectively, and transition objects by "stars". 
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Fig. 3. Comparison between equivalent widths measured in the 
FEROS solar spectrum and the SOFIN and UVES solar spectra. 
The average differences are (EWferos - EWsofin) = 1.02 ± 
2.1mA (251 lines in common) and (EWferos - EWuves) - 
0.58 + 1.9 mA. (165 lines in common) 

5. Abundance analysis 

5.1. Methods and solar analysis 

5.1 .1 . Abundances from equivalent widths 

Elemental abundances for Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, 
Zn, Y, and Ba have been determined by means of equivalent 
width measurements. The method and the atomic data are the 
same as in Bensby et al. J2003I I and are extensively described 
therein, except for the Yn and Ban lines that are listed in 
Table |4] The log ^/-values for these lines were taken from 



Pitts & Newsom J1986> for Y n, and Sneden et al. J1996t and 
Gallagher f 1967^ for Ban. We have not taken hyperfine struc- 
ture or isotopic shifts into account when measuring equivalent 
widths for the Ba n lines. The structure of the Ba n lines that 
we have used are dominated by strong central peaks, contain- 
ing the even isotopes, and smaller peaks on the sides caused 



NLTE effects while Fen lines generally are not (e.g. Thevenin 
& Idiart 1999 Gratton et al. [D^ . Figure |2 shows the dif- 
ference [Fei/Fen] versus [Fe/H], Teft, \ogg, and There are 
slight indications that the Fe i abundances come out too low for 
the stars with the highest Teff and (Figs.|2j) and d). The effect 
seems to be small (< 0.05 dex) and the number of stars at these 
higher values of Teft and are too low to allow us to delve 
deeper into this. Generally there are no significant differences 
or trends with either of the atmospheric parameters, which in- 
dicates that NLTE effects for Fe i are not severe for the majority 
of our stars. This is also true for the stars analyzed in Bensby et 
al. (2003i. The derived stellar parameters are listed in Table|3l 
for the 36 new stars. The parameters for the remaining 66 stars 
are given in Bensby et al. (■2003J . 



Table 4. Atomic line data. Column 1 gives the element; col. 2 
the wavelength, col. 3 the lower excitation potential; col. 4 
the correction factor to the classical Unsold damping constant; 
col. 5 indicates if the broadening by collisions have been taken 
from Anstee et al. {199?'}, Barklem & O'Mara dl^^Tll^^St . 
and Barklem et al. (1998 2000 1 (indicated by an "S") in- 
stead of the classical Unsold broadening (indicated by an "U"). 
Column 6 gives the radiation damping constant; col. 7 gives the 
logg/-values; col. 8 gives the references for the log §/- values. 
The full table is available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 
via cdsweb.u-strasbg . fr/cats/cgi-bin/qcat?J/A+A/ 
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Fig. 4. Synthetic and observed solar (FEROS) spectra for: a) the Euii at 4129 A, and b) the Euii line at 6645 A. Five synthetic 
spectra with different Eu abundances, in steps of 0.03 dex have been plotted for each hne. The hyperfine components for the two 
Eu isotopes are indicated as well as other important lines in the regions. The bottom panels shows the differences between the 
observed and synthetic spectra. 



by the hyperfine components of the odd isotopes (Karlsson 
& Litzen' 1999V The shifts of the even isotopes in the central 
peaks are indeed small, and are not resolvable in our spectra. 
Hence the Ban lines have essentially Gaussian line profiles. 

Figure |^ shows a comparison between the equivalent 
width measurements in the FEROS solar spectrum and the 
SOFIN solar spectrum. For the 25 1 lines in common between 
the two spectra there is a slight offset present. The FEROS 
equivalent widths are on average 1.02 + 2.1 mA larger There 
is no obvious reason for this. It can, however, be due to the 
different resolutions of the two spectrographs. In the SOFIN 
spectra, with their higher resolution, it is easier to avoid small 
blends that are not detected in the FEROS spectra, and there- 
fore the equivalent width could be on average smaller The 
SOFIN spectra in general also have higher S /N ratios which 
could lead to a lower placement of the continuum, as compared 
to the more noisy spectra from FEROS, when doing the mea- 
surements. Figure|3j' shows a similar comparison between the 
FEROS and the UVES' solar spectra. The difference is smaller 
than between SOFIN and FEROS, but the trend persists, i.e. 
that the FEROS equivalent widths are slightly higher 

The derived solar elemental abundances are tabulated in 
Tablelsjfor the different spectrographs. As expected the SOFIN 
abundances are somewhat lower In order to put all observa- 
tions on a common baseline we subtract the difference between 
the solar abundances we derive and the standard photospheric 
abundances as given in Grevesse & Sauval (1998). It should 
be emphasized that this normalization is done individually for 
each set of abundances for the different spectrographs. 



' We did not obtain a solar spectrum with UVES during our 
observing run. Instead we have used the UVES solar spec- 
trum that is available on the web: www.eso.org/observing/ 
dfo/quality/UVES/pipeline/solar.spectrum.htinl 



5.1 .2. Synthesis of the forbidden oxygen line at 6300 A 

Synthesis of the forbidden oxygen line at 6300 A has been done 
for the SOFIN spectra using the same methods and the same 
atomic data as in Bensby et al. ( 2004a i, apart from one thing: 
While we in Bensby et al. ( 2004a i used elliptical lineprofiles 
to model the combined broadening of rotation and macrotur- 
bulence we have here used radial-tangential (Rad-Tan) profiles 
instead. The difference is small and is only reflected in a slight 
shift in the absolute abundances. Since we always normalize 
our derived abundances to our own solar abundance this effect 
is of less importance as is the choice of the standard abundances 
in col. 3 in Table|5] 

We do, however, note that the solar oxygen abundance that 
we derive from the forbidden [Oi] Hne at 6300A (e(0) = 8.71 
from CES spectra and e(0) - 8.74 from SOFIN spectra are 
in good agreement with the new solar photospheric values by 
Asplund et al. ( 2004 1. They derived e(0) = 8.69 using 3D mod- 
els and e(0) - 8.73 using the MARCS model. A strict compar- 
ison between our study and theirs is, however, not straightfor- 
ward since we have used a slightly lower log gf -value for the 
[O i] line and slightly higher log g/- values for the blending Ni i 
lines. While we used log gf = -9.82 for the [O i] line they used 
\og gf = -9.72. For the blending Nil lines we used the new 
laboratory value, log gf - -2.11, from Johansson et al. J2003t 
which is split into \og gf = -2.275 for the ^^Ni component 
and \og gf = -2.695 for the ^°Ni component (see Bensby et 
al. 2004a I. Asplund et al. ©04) used \og gf = -2.31 taken 
from their previous work (Allende Prieto et al. l2001> . 

5.1 .3. Synthesis of the europium lines 

Determination of europium abundances have been done by syn- 
thesis of the Eun line at 6645 A (in the FEROS spectra) and 
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Table 5. Elemental abundances from the solar analysis. The first column gives the element and degree of ionization. An asterisk 
in the second column indicates that the log ^/-values for these lines are astrophysical (originating from Bensby et al. l2003l l. and 
an asterisk within parenthesis that a part of the lines have astrophysical log ^/-values. In the lower part of the table, col. 2 is used 
to indicate the wavelength of the spectral line. The third column gives the standard solar photospheric abundance as given in 
Grevesse & Sauval J1998> . Columns 4-6 give our solar analysis based on the FEROS spectra (see Bensby et al. l2003> . columns 
7-9 our solar analysis based on the SOFIN spectra, and columns 10-12 our solar analysis based on the UVES spectra. For each 
study we give the number of lines that were analyzed (A^iines), the mean abundance from these lines (e(X)), and the difference 
(Diff.) compared to the standard photospheric value in column 3. 
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the Eun hne at 4129 A (in both FEROS and SOFIN spec- 
tra). The synthesis was done in the same manner as for the 
[Oi] line at 6300 A. Europium has two isotopic components. 
In the solar system 47.8 % of Eu is in the form of '^'Eu and 
52.2 % in the form of '^^Eu. Europium also shows large hyper- 
fine splitting which have to be taken into account in the abun- 
dance determination. Linelists and log ^/-values for the hyper- 
fine components of the Eu lines have been kindly provided by 
C. Sneden and are the same as those used in the study by Lawler 
et al. J2001> . All atomic data are given in Table |4] 

Examples of the synthesis of the two Eu ii lines are shown 
in Fig. 0] where the different isotopic hyperfine components 
are indicated. Both Eu lines are also more or less blended with 
other lines. These lines have been taken into account in the 
modelling of the spectra and are indicated in Fig.l^as well. 

The level of the continuum was in the case of the 4129 A 
line determined from the points just to the left of the Fe ii line 
at 4128.7 A and just to the right of the Fe i Hne at 4130.0 A. For 
the Eu II line at 6645 we used the continuum points at 6644.7 A 
and 6645.7 A (see Fig.|4}. For the instrumental broadening due 
to the resolution of the FEROS and SOFIN instruments we 
adopted Gaussian profiles with apropriate widths. The com- 
bined broadening due to macroturbulence and stellar rotation 



was determined by fitting radial-tangential (Rad-Tan) profiles 
to the Fell line at 4128.7 A and the Nil Hne at 6643.6 A, re- 
spectively. 

5.2. Random errors 

The effects on the derived abundances due to random (internal) 
errors were estimated in Bensby et al. (2003 1. By changing T^ff 
by +70 K, logg by +0.1, by +0.15 kms"', [Fe/H] by +0.1, 
and the coiTection term to the Unsold approximation of the Van 
der Waals damping by +50 %, the effects were studied on four 
stars. The average of the total random error from these four 
stars are listed in Table|6lfor various abundance ratios. It should 
be noted that these estimates are made under the assumptions 
that the different error sources are uncorrected, which might 
not be completely true. Errors in the effective temperature will 
for instance also show up in the plot of Fe i abundances versus 
line strength, i.e. in the tuning of the micro turbulence. Hence, 
if this erroneous Teff was used in an analysis the researcher 
would adjust to achieve equilibrium and in this way proba- 
bly partly compensate the erroneous T^ff with an erroneous ^t- 
However, the equilibrium would not be as good as if the correct 
parameters had been used. Therefore the total errors in Table|6l 
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Table 6. Estimates of the effects on the derived abundance ratios due to internal (random) errors. The estimates are the average 
of four stars (see Bensby et al. l2003> . 
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should be treated as maximum internal errors. That these really 
are maximum errors is further reflected in the tight abundance 
trends we obtain where the scatter within each stellar popula- 
tion are lower than these estimated errors. 

5.3. Systematic errors: Comparison to other studies 

Systematic errors are more diflicult to examine. In Bensby et 
al. J2003> we compared our solar equivalent widths (measured 
in a FEROS spectrum) to those in Edvardsson et al. (1993i and 
found good agreement. The comparisons between our SOFIN 
and UVES equivalent widths to our FEROS equivalent widths 
(see Fig. |5} showed that there are small offsets present. Since 
we normalize the abundances from the different spectrographs 
separately (see Sect. 15. U we do not expect there to be any off- 
sets in the derived abundances between the different data sets. 

In Bensby et al. (.2003 ) we compared our derived abun- 
dances for a few stars to abundances from other works and 
found, generally, good agreement. In the present sample (i.e., 
all 102 stars) we have four stars in common with Reddy et 
al. (l2003t : HIP 11309 (HD 15029), HIP 85007 (HD 157466), 
HIP 92270 (HD 174160), and HIP 118010 (HD 224233). In 
Fig. Is] we compare our abundances with theirs. For three of 
the stars the differences in [X/Fe], X being any of the elements 
considered, are small apart from for one or two of the elements. 
For HIP 1 1309 the mean difference is +0.005 ± 0.048 dex, for 
HIP 85007 the mean difference is +0.019 + 0.063 dex, and for 
HIP 92270 the difference is -0.016 ± 0.058 dex. We note that 
for these three stars it is essentially two elements that con- 
tribute to the scatter, Y and Ba. These two elements also show 
systematic differences between the two studies in that we al- 
ways derive larger Ba abundances and smaller Y abundances 
than Reddy et al. (2003). If Ba and Y elements are removed 
from the calculation the resulting mean differences and scat- 
ters become +0.005 + 0.039, +0.020 + 0.036, -0.024 + 0.045, 
for HIP 11309, HIP 85007, and HIP 92270, respectively. For 



HIP 118010, however, the scatter around the mean difference 
is larger, 0.091 dex, and appear to be more random in nature. 
The scatter is not decreased when Ba and Y are removed. The 
reason for this difference lays in that we use an effective tem- 
perature that is 200 K lower than the one Reddy et al. (I2003t 
use. We have adopted a T^ti of 5795 K and Reddy et al. (l2003t 
use 5609 K. Both studies use the same surface gravity, 4. 17 dex. 
This combination results in that we derive [Fe/H] = -0.07 and 
they arrive at -0.24 dex. If we use Table|6lto estimate the cor- 
rection for this difference we arrive at abundances that are very 
similar to those of Reddy et al. ( I2003t . For the other three stars 
the stellar parameters are identical, within the errors, between 
the studies. From this comparison we conclude that our data 
and Reddy et al. (I2003t are in good agreement and when dif- 
ferences occur they can be understood and corrected for. This 
means that it is possible to combine the results from our stud- 
ies with those of Reddy et al. (1200 3> when there is a need for 
large data samples, e.g. when comparing models of chemical 
evolution to data. 

Our stellar sample have eight stars in common with 
the studies by Mashonkina & Gehren (2000 2001) and 
Mashonkina et al. (2003 ). In Fig.|5]we show a comparison be- 
tween our Fe, Ba, and Eu abundances with theirs for these eight 
stars. Except for two Ba abundances (HIP 699 and HIP 107975) 
the differences are small. For [Fe/H] the difference is 0.00 ± 
0.04 dex, for [Ba/Fe] the difference is 0.10 + O.lOdex, and for 
[Eu/Fe] the diff^erence is 0.01 ± 0.04 dex. Excluding HIP 699 
and HIP 107975 the diff'erence for [Ba/Fe] shrinks to 0.05 ± 
0.02 dex. The reason for the deviating Ba abundances for these 
two stars is difficult to resolve. If we look at the stellar parame- 
ters we see that our effective temperatures and surface gravities 
are higher than what Mashonkina and collaborators have used. 
For HIP 699 we have T^ff = 6250 K and logg = 4.19 while 
they have T.^ = 6150K and logg = 4.06, and for HIP 107975 
we have Tetj - 6460 K and logg = 4.06 while they have 
Teff = 6310K and logg - 3.94. From Table|6lwe see that this 



T. Bensby et al.: a-, r-, and .s-process element trends in the Galactic thin and thick disks 



9 



0.3 



0.2 



0.1 



-0.1 



-0.2 



A HIP 11309 DHIP 85007 i^HlP 92270 •HIP 118010 



□ 
A 



■0 if- • 



A i^r # i^lr ^ 
• □ 



□ 



Fe Na Mg Al Si Ca Ti Cr Ni Zn Y Ba 



Fig. 5. Comparison of our abundances with Reddy et al. (120031 
for four stars in common. The differences for each individual 
star are marked by symbols as indicated in the figure. The dif- 
ferences are defined as [X/Fe]This Work - [X/Fe]Reddy except for 
Fe which is [Fe/H] This work - [Fe/H]Reddy Note that for four 
elements, Ca, Ti, Cr, Ni, the differences between two stars are 
identical and hence we only see three symbols. 

does not resolve the discrepancies for these stars. If we were to 
lower our effective temperatures our [Ba/Fe] ratios would ac- 
tually increase and make the differences larger A lowering of 
the surface gravities would on the other hand lower our [Ba/Fe] 
ratios but not sufficiently. A lowering of log g by 0. 1 dex would 
only result in a 0.01 dex lowering of [Ba/Fe]. It is worth noting 
that none of these stars have abundances that make them devi- 
ate from the abundance trends that the rest of our stars outline 
(see Sect. l7.7l for stars that do). The otherwise good agreement 
to the works by Mashonkina and collaborators should justify 
combinations of stars and elemental abundances from their and 
our works when needing larger data sets. From the good agree- 
ment between our and Mashonkina and collaborators' Eu abun- 
dances we estimate systematic errors to be small for Eu, around 
or below 0.05 dex in both [Eu/Fe] and [Eu/H]. 

6. Stellar ages 

In Bensby et al. (I2003> we used the Salasnich et al. (^000 1 



and Girardi et al. (I2000> isochrones to estimate the stellar 
ages. The Yoshii-Yale (Y^) isochrones (Kim et al 120021 Yi 
et al. 1200 U are distributed together with interpolation rou- 
tines that makes it possible to construct a set of isochrones 
with any metallicity and a-enhancement. Other isochrone sets 
are tabulated for fixed values of these parameters. These fixed 
values may not necessarily coincide with the analyzed data. 
Since a-enhancement varies for our stars we have opted to 
use the Y^-isochrones. At sub-solar [Fe/H] we used different 




-0.2 



<V <\> rs° (W' rf^ (h 

^ ^ # # # # ^ ^ 



Fig. 6. Comparison of abundances for eight stars that 
we have in common with the studies by Mashonkina & 
Gehren (2000 2001) and Mashonkina et al. 



( I2003> . The dif- 
ferences are defined as our abundances minus theirs. Note that 
the Mashonkina studies have no Eu abundances for two of the 
stars: HIP 93185 and HIP 107975. 

Table 7. Metallicities and a-enhancements for the Y^ 
isochrones that were used in the age determination. 



[Fe/H] 



[a/Fe] 
Thin Thick 



[Fe/H] 



Thin Thick 



-0.90 
-0.80 
-0.70 
-0.60 
-0.50 
-0.40 
-0.30 



+0.15 
+0.13 
+0.12 
+0.10 
+0.07 



+0.30 
+0.30 
+0.30 
+0.30 
+0.30 
+0.30 
+0.20 



-0.20 
-0.10 
0.00 
+0.10 
+0.20 
+0.30 
+0.40 



+0.05 
+0.03 
+0.03 
+0.03 
+0.03 
+0.03 
+0.03 



+0.15 
+0.10 
+0.03 
+0.03 
+0.03 
+0.03 
+0.03 



ff-enhancements for thin and thick disk stars according to our 
observations (see Table 0. The most likely age was then es- 
timated for each star from isochrones plotted in the Teff - My 
plane, using Hipparcos parallaxes and our spectroscopic tem- 
peratures. Lower and upper age limits were estimated from the 
plots as well by taking the errors in the parallaxes and effec- 
tive temperaures into account. Stellar ages were determined for 
the new stars in this study and, in order to get a consistent age 
determination for the whole stellar sample, those in Bensby et 
al. (l2003i as well. The ages and their lower and upper limits 
are given in Table|F] Other methods to derive stellar ages from 
isochrones exist. However, our age estimates are virtually iden- 
tical to those obtained with more sophisticated methods. Our 
method, on the other hand, probably underestimates the un- 
certainties of the derived ages (Rosenkilde J0rgensen, private 
comm.). 

The mean ages of the thin and thick disk samples (includ- 
ing the new age determinations for the stars from Bensby et 
al. lStKBt are 4.3 + 2.6 Gyr and 9.7 + 3.1 Gyr, respectively. The 
mean age for the four stars with intermediate kinematics is 
6.7 ±2.0 Gyr. 
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Table 8. Age estimates for the stars in this study and 
those in Bensby et al. (2003 1. The minimum (Min age) 
and maximum (Max age) ages are based on the uncer- 
tainties in the Hipparcos parallaxes and +70 K in the ef- 
fective temperatures. The full table is available via anony- 
mous ftp at cdsarc.u-strasbg.fr (130.79.128.5) or via 
cdsweb.u-strasbg . £r/cgi-bin/qcat?J/A+A/ 



Table 9. Derived abundances relative to hydrogen. Each 
element has three columns, mean abundance ([X/H]), stan- 
dard deviation of the mean abundance (ct^x/h]), and the 
number of spectral lines (A^x) that were used in the abun- 
dance analysis (Cols. 4 and onwards). The abundances have 
been normalized with respect to the Sun (see Sect. 15.11 and 
Table |5|i. Column 1 gives the Hipparcos number; col. 2 
indicates with which spectrograph the star was observed 
(S=SOFlN, F=FEROS, U=UVES); col. 3 indicates if the 
star belongs to the thin disk (Mem.=l), the thick disk 
(Mem. =2), or a "transition object" (Mem. =3). The full ta- 
ble is available in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
cdsweb .u-strasbg . fr/cgi-bin/qcat?J/A+A/ 
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Figure0shows [Fe/H] as a function of age for all 102 stars. 
For the thick disk stars there might exist a relation between 
age and metallicity. Thick disk stars with [Fe/H] < -0.4 have 
a mean age of 11.6 + 2.9 Gyr, those with -0.4 < [Fe/H] < 
-0.2 have a mean age of 8.1 ± 2.7 Gyr, and those with [Fe/H] > 
-0.2 have a mean age of 7.7 + 1.4 Gyr This indicates that star 
formation could have continued in the thick disk for quite some 
time, up to about 2-3 Gyr This conclusion is however uncertain 
due to the large spread in the stellar ages and especially to the 
rather small stellar sample we have here. The potential trend 
between age and metallicity in the thick disk is, however, very 
similar to the results we find in a study of ages and metallicities 
(based on Stromgren uvby photometry) of a larger sample of 
thick disk stars (Bensby et al. i2004bi) . In that study we found a 
possible age-metallicity relation in the thick disk and also that 
it might have taken 2-3 Gyr for the thick disk to reach [Fe/H] « 
-0.4 and a further 2-3 Gyr to reach solar metallicites. 

We note that at the highest metallicities ([Fe/H] > 0.2) there 
is a lack of young stars (ages lower than 3 Gyr). Given that there 
is ongoing star formation in the metal-rich thin disk today this 
is clearly not a representative picture. This apparent trend is 




5 10 15 20 



Age [Gyr] 

Fig. 7. [Fe/H] versus age for the 36 stars from this study and 
the 66 stars from Bensby et al. (2003). Thin disk and thick 
disk stars are marked by open and filled circles, respectively. 
Transition objects are marked by asterisks. 



probably due to selection effects in our sample, i.e. only con- 
tructed from F and G dwarf stars and not including earlier type 
stars, and not a feature of the Galactic thin disk. 



7. Abundance results 

7.1. Strengthening the a- and iron-peak element 
trends 

In Fig. Owe show the abundance trends relative to Fe for all 
102 stars. The new stars from the northern sample confirm and 
extend the trends that we presented in Bensby et al. J2003I 
'2004a'l. No major novelties are found so these trends will only 
be briefly described and the reader is directed to Bensby et 
al. (120031 12004a> for further discussions and comparisons to 
other works. 

O, Mg, Si, Ca, and Ti: Thin and thick disk stars are clearly 
separated and show tight and distinct trends. The down-turns 
(or "knees") in the [O/Fe] and [Mg/Fe] trends for the thick disk 
stars are especially prominent (see Fig.|8^ and c). There is no 
doubt that the locations of these "down-turns" are at [Fe/H] =s 
-0.4 after which [O/Fe] and [Mg/Fe] decline toward solar val- 
ues. At lower [Fe/H] the trends are essentially flat showing con- 
stant values of [Mg/Fe] ~ +0.35 and [O/Fe] ~ +0.4. 

Although not as prominent as for O and Mg, these features 
are also clearly present for the other three a-elements; Si, Ca, 
and Ti. 



HIP Mem. Age Min age Max age 



699 
910 



3.8 
5.5 



3.4 
5.0 



4.4 
5.8 
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Fig. 8. Elemental abundances relative to Fe. Dotted lines indicate solar values. Thin disk and thick disk stars are marked by 
empty and filled circles, respectively. Stars from Bensby et al. (.2003. ■2004aJ are marked by circles and stars from the new 
northern sample by triangles. Transition objects are marked by asterisks. 



Na and Al: The appearance of the thin and thick disk [Na/Fe] 
trends (Fig. (Sji) show an inverted behaviour compared to the 
a-elements. Instead, the stars in thin disk seem to be more 
abundant in Na than those in the thick disk. At first glance 



this [Na/Fe] trend appear to be different from that in Bensby et 
al. 12003 ) where we found that [Na/Fe] trends for the thin and 
thick disk appear to be merged. The tighter [Na/Fe] trends in 
this study (more stars are now used to trace the trends) indicate 
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that also the Na abundances are distinct between the two disk 
populations, even if not as well separated as for the a-elements. 

The Al trends show the same type of trends as the a- 
elements (see Fig. |8ji). This supports our findings in Bensby 



et al. J2003> that Al and the or-elements are produced in the 
same environments and have been dispersed into the interstel- 
lar medium on the same time-scales, i.e., the SN 11 events. Also, 
Mc William ( 1997 1 noted that Al, from a phenomological point 
of view, can be classified as an or-element. 

Cr and Ni: Cr varies tightly in lock-step with Fe (see Fig.|8h). 
No trends are seen and the thin and thick disk stars are well 
mixed which strongly emphasize the common origin for Cr and 
Fe. 

In Bensby et al. i2003l l we found that the Ni and Fe below 
solar metallicities evolve roughly in lockstep (i.e., [Ni/Fe] x 0). 
At [Fe/H] ~ the [Ni/Fe] trend then showed a prominent up- 
turn that had not been seen in previous studies. This deviation 
from a flat [Ni/Fe] trend has an impact on the oxygen abun- 
dances that are derived from the forbidden [O i] line at 6300 A 
since this hne is heavily blended by two Ni i lines (see Bensby 
et al. l2004aV In the new sample there are only two stars with 
[Fe/H] > 0. They do, however, follow the same trend as found 
in Bensby et al. (2003 1 and show an increased [Ni/Fe] (see 
Fig. |8l)- The now larger number of stars at [Fe/H] < also 
indicates that it is possible that the [Ni/Fe] trend at these metal- 
licities actually is not flat. We see a slight overall decrease in 
[Ni/Fe] when going to higher [Fe/H], and at [Fe/H] = there is 
an underabundance of Ni relative Fe of about 0.05 dex. There 
is also a weak tendency that the thick disk stars are more abun- 
dant in Ni than the thin disk stars. 

Zn: The [Zn/Fe] trend is tight and in accordance with Bensby 
et al. J2003> (see Fig.|8]). This is somewhat surprising since the 
trend for the new stars is based on one Zn i line only. Further, 
this is also the hne that we rejected from further analysis in 
Bensby et al. ( I2003> since it was suspected to have a hidden 
blend that was growing with metallicity. For metallicities below 
[Fe/H] X! we, however, found that the blend should have less 
influence. This is most likely the reason for the good agreement 
between our [Zn/Fe] trends. 

Comparing our thin disk [Zn/Fe] trend with Reddy et 
al. (20031 we see that in the range -0.4 < [Fe/H] < -0.2 
their stars have [Zn/Fe] in the range -0.1 dex to -1-0.2 dex. This 
is higher than what we see for our thin disk stars that have 
[Zn/Fe] in the range -0.1 dex to Odex in the same metallic- 
ity bin. Our thick disk stars have [Zn/Fe] in the range Odex 
to +0.2 dex which means that by combining the [Zn/Fe] trends 
for our thin and thick disks we would see the same spread in 
[Zn/Fe] as Reddy et al. (2003 1. However, as we saw in Sect. l5.3l 
there seem to be an offset of about 0.05 — 0.10 dex between our 
[Zn/Fe] abundances and those in Reddy et al. (2003 1. Taking 
this into account will put our thin disk [Zn/Fe] trend on the 
same level as the one in Reddy et al. ( 2003), or vice versa. But, 
it will not give any insight into why the thin disk [Zn/Fe] trend 
in Reddy et al. (t2003j show a larger scatter than what we see 
in our [Zn/Fe] trend. It is probably due to the analysis in which 
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Fig. 9. [O/Fe] versus Zmax for the 36 stars from this study and 
the 66 stars from Bensby et al. ('2003'). Thin disk and thick 
disk stars are marked by empty and filled circles, respectively. 
Transition objects are marked by asterisks. 
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Fig. 10. Abundance trends in the thick disk for oxygen. The 
thick disk sample has been divided into two sub-groups: stars 
that have Z^ax > 500 pc (open squares); and stars that have 
Zmax < 500 pc (open stars). 



only one or two spectral lines are analyzed, which unevitably 
leads to larger internal errors unless extreme care is taken. 

7.2. [a/Fe] at different Z^ax 

Since our thick disk stellar sample is far from complete and is 
biased towards higher metallicities we can not use it to probe 
for vertical gradients in the thick disk metallicity distribution. 
However, it can be used to investigate if there are differences 
in the abundance trends at various heights above the Galactic 
plane. If the trends are similar this would indicate that the thick 
disk stars come from a stellar population that initially was ho- 
mogeneous and well mixed. 
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The maximum vertical distance (Z^ax) ^ star can reach 
above the plane can be estimated from (L. Lindegren 2003, 
priv. comm.): 

Zmax = + (k ■ Wi^sr)\ (1) 

where 

k = (85 Myr) / (2 ■ tt) (pc / Myr) / (km / s), (2) 

and Z is the present distance of the star from the Galactic 
plane (in pc) and Wlsr is its present velocity in the Z-direction. 
Equation Q assumes that stars in the solar neighbourhood os- 
cillates harmonically in the Z-direction with a period of 85 Myr, 
independent of their other motions. Since the density of stars 
decreases with |Z|, Eq. Q will underestimate Z„iax for stars that 
have large Wlsr velocities. The relationship is however suffi- 
cient for our purposes. 

Figure|9]shows the [O/Fe] ratio as a function of Z^ax for our 
sample. The thin disk stars are all confined to within ~ 300 pc 
of the Galactic plane, whereas the thick disk stars move in or- 
bits reaching vertical distances up to 1 kpc or more. Dividing 
the thick disk sample at ~500pc give two sub-samples of ap- 
proximately equal sizes. 

As can be seen in Fig.^Jthe [O/Fe] trends for the two sub- 
samples are the same. An important point is that the "knee" is 
located at the same [Fe/H] for both sub-samples. If the Galactic 
thick disk formed in a fast dissipational collapse, with a pro- 
posed time scale of ~ 400 Myr (see, e.g., Burkert et al. ll992l l. it 
is likely that the position of the "knee" would differ in the two 
sub-samples. SNll that have a time-scale of typically 10 Myr 
would then have time to enrich the interstellar medium with 
even more of the a-elements at higher [Fe/H] in the thick disk 
sub-sample that formed closer to the Galactic plane. The in- 
variance of our abundance trends with distance from the plane 
instead indicates that the thick disk stellar population was well 
mixed before it got kinematically heated. This can for example 
be accomplished if the stars in a pre-existing old thin disk got 
kinematically heated by the tidal interaction with a companion 
galaxy that either merged with, or passed close by, the Galaxy. 

7.3. The r- and s-process element abundance trends 

Figure II Ih-c shows our results for the r- and .s-process ele- 
ments Y, Ba, and Eu with Fe as the reference element, and 
Fig. II IH-f with oxygen as the reference element. 

Eu: Eu is supposed to be an almost pure r-process element 
(approximately 94 % r-process and 6 % i-process, Arlandini et 
al.D999i which means that it mainly should come from SN 11 
where the neutron flux is sufficient for the r-process to occur 
(i.e. the neutron density is so high that the neutron-capture 
timescale is much smaller than the /?-decay timescale). The 
[Eu/Fe] trend in Fig.lllh is very similar to the [O/Fe] trend in 
Fig- lit- The thick disk [Eu/Fe] trend also shows a turn-over at 
[Fe/H] X! -0.4, and the thin disk [Eu/Fe] trend shows the same 
shallow decline over the whole range in [Fe/H]. The downward 
trend that we found for [O/Fe] at [Fe/H] > is also present in 
[Eu/Fe] but with a larger scatter. The generally good agreement 



between Eu and oxygen, which is further illustrated in Fig. II If 
where we plot [Eu/O] versus [O/H], indicates that these two 
elements indeed originate from the same type of environments. 

Our [Eu/Fe] trend for the thin disk is in good agree- 
ment with previous studies (Woolf et al. 119951 Koch & 
Edvardsson 120021 Mashonkina & Gehren '200 1^, and our 
thick disk [Eu/Fe] trend is in agreement with Mashonkina & 
Gehren (1200 U for the metalUcities where their thick disk stars 
overlap with ours (i.e. [Fe/H] < -0.3). The continuing decline 
that we see in [Eu/Fe] for the thick disk for [Fe/H] > -0.3 is, 
on the other hand, new. 

Ba and Y: The s-process contributions to the solar composi- 
tion is for Ba 8 1 % (Arlandini et al. 1999 but see also Travaglio 
et al. 1999'1 and for Y 74% (Travagho et al-EOOS- In the s- 
process the neutron flux is low, which means that the radioac- 
tive isotopes will have time to /?-decay between the neutron- 
captures. Probable sites for the s-process are the atmospheres 
of stars on the asymptotic giant branch (AGB stars) (e.g., 
Busso et al. 1999 1. If this is the case the enrichment of the s- 
process elements to the interstellar medium probably occur on 
a timescale similar to elements originating in SN la. 

The [Y/Fe] and [Ba/Fe] trends are different for the thin and 
thick disks (Fig.ll Ih and b). Especially [Ba/Fe] is distinct and 
well separated for the two disks. For the thick disk stars the 
[Ba/Fe] trend (Fig.ll lb) is flat, lying on a solar ratio. [Y/Fe] for 
the thick disk shows a larger scatter and has a flat appearance 
with underabundances between and -0.2 dex (Fig.ll Ih). The 
thin disk [Ba/Fe] trend shows a prominent rise from the lowest 
[Fe/H] until reaching solar metallicities, after which it starts to 
decline. The [Y/Fe] trend for the thin disk is similar but shows 
a considerably larger scatter and not such a well-defined trend 
as that for [Ba/Fe]. 

7.4. Evolution of Ba and Eu 

In Fig.[21we compare our Eu and Ba data with predictions for 
pure r-process composition ([Eu/Ba] = 0.7, calculated from 
the yields given in Arlandini et al. 1 19991 1 and a solar mixture 
of r- and s-process contributions (i.e., [Eu/Ba] = 0). Also in- 
cluded in this plot are low-metallicity halo stars (giants) from 
Burris et al. <2000l l. The thin disk shows a solar system mix for 
all metallicities while the thick disk has not yet experienced the 
full contribution of s-processed material from low mass AGB 
stars, i.e., it is closer to the pure r-process line. 

A first tentative interpretation of these results is that star 
formation went on long enough in the thick disk so that AGB 
stars started to contribute to the chemical enrichment, but only 
just long enough that a solar system mix was reached. After the 
formation of stars stopped in the thick disk the remaining gas 
settled into a new thinner disk. Most likely, fresh material of 
lower metallicity was accreted before star formation started in 
what is today's thin disk. The relative r- and s-process contri- 
butions will not change by this dilution if the infalling material 
is pristine, so in this case the first thin disk stars to form will 
retain the mixture that was at the end of star formation in the 
thick disk. The absolute abundances of Ba and Eu in the thin 
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Fig. 11. Elemental abundances relative to Fe. Dotted lines indicate solar values. Thin disk and thick disk stars are marked by 
empty and filled circles, respectively. Transition objects are marked by asterisks. 



disk will, however, be shifted towards lower values. If on the 
other hand, the infalling material has experienced enrichment 
of r- and/or s-process elements the mixture should change. Our 
data appear to indicate that this has not been the case and hence 
that the infalling material was most likely primordial. 

This is further illustrated in Fig.[0]where we plot [Eu/Ba] 
versus [Fe/H] for our stellar sample only. The most metal-rich 
thick disk stars (at [Fe/H] x 0) and the most metal-poor thin 
disk stars (at [Fe/H] ^ -0.7) have approximately the same 
[Eu/Ba] ratio. So while pristine material falls into the disk the 
gas (from which the thin disk stars form) gets more metal-poor, 
the [Eu/Ba ratio is preserved. 

7.5. Nucleosynthesis in low mass AGB stars 

An interesting result is found when studying the trend of 
[Y/Ba] versus [O/H], see Fig.[l4| For both the thin and thick 
disks this trend is first flat but after solar metallicity a gentle 
upward trend is seen. This could be explained as a metallic- 
ity effect in AGB nucleosynthesis (Busso et al..2001 ). Figure 1 
in Busso et al. (2001) shows how the relative production of 
light (e.g. Y) and heavy (e.g. Ba) i-process elements change 
as a function of metallicity. For metallicities below solar the 
[Y/Ba] ratio is roughly flat, i.e. Y and Ba are produced in the 



same ratio in the low mass AGB stars. However, around, or 
slightly above, solar metallicity this balance changes such that 
the lighter i-process elements are favoured over the heavy s- 
process elements. Hence [Y/Ba] increases. 

7.6. Transition objects 

In our stellar sample we have four stars whose kinematical 
properties lie in between the definitions of the thin and thick 
disk populations that we have used (see Table Depending 
on the choice of the normalization of the thick disk density 
in the solar neighbourhood they will alter their classifications 
as either thin disk or thick disk stars (see discussion in the 
Appendix). When looking at their chemical compositions we 
see that they not only are intermediate in terms of kinematics 
but also in terms of abundances at a given [Fe/H] (see Figs.|8l 
and ll li . Which stellar population these stars should belong to, 
or if they form a distinct population by their own, can only be 
investigated with a larger sample of such stars. 

7.7. Deviating stars - outliers 

From our abundance plots it is evident that not all stars fol- 
low the trends as outlined by the majority of stars in the two 
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Fig. 12. An illustration of the relative contributions from the r- 
and i-process to the elements Eu and Ba. The dotted line shows 
the pure r-process contribution to Eu and Ba while the full line 
shows the solar system mix of r- and i-process contributions 
(Arlandini et al. I1999> . Thin and thick disk stars are marked 
by open and filled circles, respectively. Transition objects are 
marked by asterisks (*), and 'open stars' data from Burris et 
al. (120001. 



CQ 



3 




-0.4 -0.2 
[Fe/H] 

Fig. 13. [Eu/Ba] versus [Fe/H]. The dashed lines show the pure 
r-process ratio which is [Eu/Ba] =0.7 and has been calculated 
from the yields given in Arlandini et al. ( il999< . 

disk populations. There are, in particular, three stars (HIP 2235, 
HIP 15510, and HIP 16788) that seem to show suspiciously 
high abundances in either or both of Y and Ba (see Figs. II Ih. 
b, d, and e). To enable a direct comparison we list their chemi- 
cal properties in Table 1 101 and discuss each of them in turn. 
HIP 2235 is a thin disk star and has spectral type F6V (accord- 
ing to the Simbad database) and shows high over-abundances 
in both Ba and Y that is not recognized in any of the other thin 
disk stars (or the thick disk stars). This s-enhancement can be 
due to that s-enriched material has been transferred from from 



Table 10. Abundances for deviating stars. Each star has three 
columns; abundance; line-to-line scatter (Icr standard devia- 
tion); and (in parenthesis) number of lines that were used to 
derive the abundance. 





HIP 2235 


HIP 15510 


HIP 16788 




Thin disk 


Thick disk 


Thick disk 


[Fe/H] 


-0.28 0.07 (65) 


-0.41 0.08 (69) 


-0.32 0.07 (78) 


[O/H] 




-0.01 0.00(1) 


-0.27 0.00(1) 


[Eu/H] 


-0.19 0.00(1) 




-0.06 0.00(1) 


[Ba/H] 


0.54 0.06(4) 


-0.56 0.08 (3) 


0.63 0.06(4) 


[Y/H] 


0.31 0.14(3) 


0.09 0.00(1) 


0.45 0.12(4) 



a companion star into the stellar atmosphere. The companion 
star could have been an AGB star, in which these elements are 



believed to be synthesized (see, e.g., Abia et al. l2002t . that now 
is invisible since it has evolved into a white dwarf. 
HIP 15510 is a G8V thick disk star and show an abnormal en- 
hancement in Y but not in Ba. It is not unlikely that this star not 
has been subject to the same type of mass-transfer as HIP 2235 
maybe has been. It is namely well possible to have s-stars with 
Y-enhancement and no Ba-enhancement and vice versa (see, 
e.g., review by Busso et al. 2004 1. However, the Y abundance 
for HIP 15510 is based on one spectral line only, making it 
highly uncertain. 

HIP 16788 is a GO thick disk star that has no luminosity clas- 
sification in the Simbad database. However, from our derived 
log g = 4.24 it is most likely also a main sequence star. It is 
highly enhanced in both Y and Ba as in the case for HIP 2235. 

Stars from the "old sample" that show deviating a- 



abundances were discussed in Sect. 9.4 in Bensby et al. J2003> 
to which the reader is referred since none of the new stars 
showed deviating abundances for these elements. 



16 



T. Bensby et al.: a-, r-, and .s-process element trends in the Galactic thin and thick disks 



8. Discussion and summary 

In this study we have presented a differential abundance analy- 
sis between the Galactic thin and thick disks for 14 elements 
(O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Zn, Y, Ba, and Eu) 
for a total of 102 nearby F and G dwarf stars (including the 
stars from our previous studies Bensby et al. 120031 |5004ai. 
The results from the 36 stars in the new sample further con- 
firms, strengthens, and extends the results presented in Bensby 
et al. (I2003ll2004a> . Results that are new in this study are those 
for the r- and s-process elements Y, Ba, and Eu, where we find 
the thin and thick disks abundance trends to be distinct and 
well defined. We also see indications of a metallicity effect in 
the AGB nucleosynthesis of Y and Ba, such that Y is favoured 
over Ba at higher [Fe/H] . Our results for Eu show that Eu abun- 
dances follow the oxygen abundances very well. This confirms 
that Eu is an element that mainly is produced in SN II. 

In our studies we have included thick disk stars with 
[Fe/H] > -0.35, which no other study have. At these higher 
metallicities we find that the [or/Fe] trends, at [Fe/H] ^ -0.4, 
turns over and decline towards solar values where they merge 
with the thin disk [a/Fe] trends. The observed down-turn (or 
"knee") in the thick disk [(r/Fe] trends at [Fe/H] ^ -0.4 can 
be interpreted as a signature of the contribution from SN la to 
the chemical enrichment of the stellar population under study. 
Massive stars (M > 10 Mq) explode as core-collapse super- 
novae type II (SN II) and enrich the interstellar medium with 
a-elements and lesser amounts of heavier elements such as 
the iron peak elements (e.g. Tsujimoto et al. 11995 1' Woosley 
& Weaver [l995> . Due to the short lifetimes of these massive 
stars they enrich the interstellar medium in the early phases 
of the chemical evolution and produce high [c/Fe] ratios at 
the lower metallicities. SNIa disperse large amounts of iron- 
peak elements into the interstellar medium and none or little 
of a-elements. Since their low-mass progenitors are expected 
to have much longer lifetimes than the SN II progenitors (e.g. 
Livio 1200 It there will be a delay in the production of Fe as 
compared to the ff-elements. Hence, when SNIa start to con- 
tribute to the enrichment, the [a/Fe] ratios will decrease. 

The fact that we see the signatures from SN la in the thick 
disk thus means that star formation must have continued in the 
thick disk for a time that was at least as long as the time-scale 
for SN la. The time-scale for a single SN la is very uncertain 
(see e.g. Livio 2001 1. However, we have seen from a study of 
ages and metallicities in the thick disk that it has taken about 
2-3 billion years for the thick disk stellar population to reach 
a metallicity of [Fe/H] = -0.4 (Bensby et al. 2004b i. Thus we 
can tentatively conclude that the SN la rate peaked at ~ 2-3 bil- 
lion years from the start of the star formation in the population 
that we today associate with the thick disk. The age-metallicity 
relation in the thick disk that we find in that study also indi- 
cates that star formation might have continued for 2-3 billion 
years after the peak in the SN la rate in order to reach roughly 
solar metallicities. The most important conclusion from this is 
that the thick disk most probably formed during an epoch span- 
ning several (> 2-3) billion years. Data from solar neighbour- 
hood stars have also shown the SN la time-scale to be as long 
as ~ 1.5 Gyr (Yoshii et al. [T996t . 



We are able to draw further conclusions about the origin 
and chemical evolution of the thick disk. The observational 
constraints for a formation scenario of the thick disk are: 

1 . Distinct, smooth, and separated abundance trends between 
the thin and thick disks. 

2. At a given [Fe/H] below solar metallicities the thick disk 
stars are more enhanced in their a-element abundances than 
the thin disk 

3. The thick disk or-element trends show signatures of enrich- 
ment from SN la. 

4. The abundance trends are identical for thick disk stars that 
reach different heights above the Galactic plane (> 500 pc 
and < 500 pc). 

5. The thick disk stars have an older mean age than the thin 
disk stars 

6. AGB stars have contributed to the chemical enrichment of 
the thick disk but not as much as to the enrichment of the 
thin disk 

7. A less well established constraint from this study is that 
there might be a possible age-metallicity relation in the 
thick disk indicating that star formation might have con- 
tinued for several billion years. 

Taking these constraints into consideration we argue that 
the currently most probable formation scenario for the thick 
disk is an ancient merger event between the Milky Way and 
a companion galaxy. In this event the stellar population of 
the thin disk that was present at that time got kinematically 
heated to the velocity distributions and dispersions that we see 
in today's thick disk. We note that recent models of hierar- 
chical galaxy formation might be able to succesfully repro- 
duce thick disks in Milky Way like galaxies and the abundance 
trends might be fully explainable also in these models (Abadi 
et al. 2003 1. 

How can we explain the trends observed in the thin disk? 
The thin disk stars, on average, are younger than the thick disk 
stars. However, the low-metallicity tail in the metallicity distri- 
bution of the thin disk stars overlap with the metallicity distri- 
bution of the thick disk stars. A possible scenario would be that 
once star formation in the thick disk stops, there is a pause in 
the star formation. During this time in-falling fresh gas accu- 
mulates in the Galactic plane, forming a new thin disk. Also, 
if there is any remaining gas from the thick disk it will settle 
down onto the new disk. Once enough material is collected, star 
formation is restarted in the new thin disk. The gas, though, has 
been diluted by the metal-poor in-falling gas. This means that 
the first stars to form in the thin disk will have lower metallici- 
ties than the last stars that formed in the thick disk. 
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Fig. A.l. CMD for the full stellar sample (12 634 stars). 

Appendix A: Kinematical criteria for selecting 
thick disl< stars in the solar neighbourhood 

When selecting thin and thick disk stars we assume that the 
Galactic space velocities (Ulsr, Vlsr, and Wlsr) for the thin 
disk, thick disk, and stellar halo have Gaussian distributions. 
The space velocities Ulsr, Vlsr, and Wlsr were calculated 
using our measured radial velocities and positions, proper 
motions, and parallaxes from the Hipparcos catalogue (see 
Eqs. A.1-A.4 in Bensby et al. 2003 1. For each star we then 
calculate the probabilities that it belong to either the thin disk 
(D), thick disk (TD), or the halo (H). By also take the fraction 
of thick disk stars in the solar neighbourhood into account, the 
final relationship for calculating the individual probabilities are 
(see also Bensby et al. 2003): 



P ^X k exp 



^LSR (Vlsr - Vasym)^ 



24 



"^"^LSR 



where P is either D, TD or H; and 
1 



(27r)3/2 



CTu cry CTw 



(A.l) 



(A.2) 



normalizes the expression; cru, cry, (Tw are the characteristic 
velocity dispersions; Vasym is the asymmetric drift; and X is the 
observed fraction of stars in the solar neighbourhood for each 
population. We calculated for each star the "relative probabili- 
ties" TD/D and TD/H. Thin disk stars were selected as those 
with TD/D < 0. 1 (i.e. at least ten times more probable of being 
a thin disk stars than a thick disk star) and thick disk stars as 
those with TD/D > 10 (i.e. at least ten times more probable of 
being a thick disk stars than a thin disk star). TD/H > 1 was 
required for both the thin and thick disk stars. 

Of the parameters that has the largest influence on the de- 
rived probability ratios, the normalization factor is the one that 



Table A.l. Characteristic velocity dispersions (cru, cry, and 
cTw) in the thin disk, thick disk, and stellar halo, used in 



Eq. ( lA.U . Vasym is the asymmetric drift. 



0"U 0"V CTW K.sym 

[kms-'] 



Thin disk (D) 
Thick disk (TD) 
Halo (H) 



35 
67 
160 



20 
38 
90 



16 
35 
90 



-15 
-46 
-220 



Table A.2. The number of stars in given TD/D intervals for 
different values on the local density of thick disk stars (Xju). 
The second column indicates the factor by which the TD/D 
ratios change when varying the normalization (with the 10% 
density as base). The corresponding CM-diagrams can be seen 
inFig.ESl 



<0.1 



* Stars 

0.1-1 1-10 



> 10 



2% 


0.2 


11781 


362 


166 


261 


6% 


0.6 


11305 


689 


238 


347 


10% 


1 


10969 


946 


285 


383 


14% 


1.5 


10623 


1215 


327 


421 



is the least well constrained. In Bensby et al. ("2003^ we used a 
value of 6 %. The lowest value 2 % was found by Gilmore & 
Reid ( 1983 1 and Chen ( 1997 1. Intermediate values ~6% were 
found by Robin et al. (i 19961 and Buser et al. ( 1999 1, and higher 
values ~ 15 % by Chen et al. (jSOOl I and Soubiran et al. ( 2002). 

We will here show that an increase of the local normaliza- 
tion of the thick disk stars from 6 % to 10 % (and consequently 
a lowering of the thin disk density from 94 % to 90 %) is moti- 
vated. Changing to 10 % will not influence the thick disk sam- 
ple in Bensby et al. (2003^ since the only effect is to raise the 
TD/D ratios by a factor of ~ 1.7 (see column 2 in Table lA.2t . 
Instead those thick disk stars will have their classicications 
strengthened. The thin disk sample in Bensby et al. (i2003t l will 
not change either since all those stars had their TD/D ratios 
well below 0.1 (typically 0.01). 

We select our stars from the same data set as in Feltzing 
et al. 12001.) and Feltzing & Holmberg (2000.). In brief this 



includes all stars in the Hipparcos catalogue (ES A [19971 that 
have relative errors in their parallaxes less than 25 % and that 



have published radial velocities (see e.g. Bensby et al. 120031 
Feltzing & Holmberg 120001 ). This sample consists of 12634 
stars. Note that known binaries have been excluded from the 
data (see Feltzing et al. 2001 ). 

The colour-magnitude diagram (CMD) for all 12634 stars 
is shown in Fig. lA.ll Figure [A2l show the CMD for three dif- 
ferent TD/D intervals: 0.1 < TD/D < 1 (i.e "low probabihty" 
thin disk stars); 1 < TD/D < 10 (i.e "low probability" thick 
disk stars); and TD/D > 10 (i.e "high probability" thick disk 
stars), with four different values of the thick disk normalization 
in the solar neighbourhood (2 %, 6 %, 10 %, 14 %). The CMD's 
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Fig. A.2. CMDs for different interval of TD/D when changing the local density of thick disk stars (2 %, 6 %, 10 %, and 14 %). 
The number of stars in each CMD is given in Table lA.21 



T. Bensby et al.: a-, r-, and .s-process element trends in the Galactic thin and thick disks 



19 



for the "high probability" thin disk (i.e. TD/D < 0.1). They are 
all essentially all hke the CMD in Fig. lA.ll 

Thick disk {TD/D > 10): The first thing that is evident from 
Fig. IA.2I is that the proposed thick disk population have a 
prominent turn-off and that the CMD essentially do not change 
when going from the lowest to the highest value for the normal- 
ization. All stars that move into this TD/D range when raising 
the normalization come from the 1 < TD/D < 10 bin. Hence 
they are always classified as thick disk stars regardless of the 
normalization. The prominent turn-off that is seen is typical for 
an old stellar population which is in concordance with the cur- 
rent beliefs of the Galactic thick disk (see e.g. Fuhrmann il998,) . 

Tliick disk (1 < TD/D < 10): The CMD's for these stars re- 
semble those for the stars with TD/D > 10. The number of 
stars in this TD/D range almost double when going to the high- 
est thick disk normalization. As in the case for the TD/D > 10 
stars not much happens at first glance when raising the thick 
disk normalization. However, when going from the 10 % to the 
14 % normalization a few stars with My < 5 and B - V < 0.5 
(i.e. in the upper left area of the CMD) populate the CMD for 
the higher normalization. These stars are most likely younger 
objects that should be attributed to the thin disk. A 14 % nor- 
malization is therefore possibly too high given the velocity dis- 
persions and rotational lags that we use (see Table lA.ll . 

Tliin disk (0.1 < TD/D < 1): A young population such as 
the thin disk should include young stars, especially the more 
massive stars that are still located on the main sequence. For 
the 2 % normalization the CMD look suspiciously like the 
CMDs for the samples with higher TD/D ratios (i.e. the thick 
disk stars). The upper left hand area is poorly populated and 
the CMD probably mainly consists of stars with "hot" thick 
disk kinematics that have a too low TD/D ratio due to a 
too low normalization. The picture is improved for the 6 % 
normalization where a substantial number of young objects 
start to populate the CMD. 

From this simple investigation we conclude that it is likely 
that a normalization of 2 % is too low and a 14 % normalization 
probably is too high for the thick disk in the solar neighbour- 
hood. Somewhere in between there is a dividing line where ob- 
viously young objects starts to populate the thick disk CMD. 
The exact value for this normalization is of course also depen- 
dent on the assumed velocity dispersions in the disks. With our 
aaumptions it is however probably located closer to 10 %. We 
have therefore used Xj£, = 10 % in the calculation of our prob- 
abilities. 

Selecting stars for abundance analysis: Where to put the 
limit where to assign a star as a "true" thick disk star is how- 
ever difficult. The most safe way is to make it independent of 
the thick disk normalization. When doing this we will disre- 
gard the 2 % normalization since it is obviously too low, and 
only consider values on X-yd between 6 and 14 %. When calcu- 



lating the TD/D ratios with the 10 % normalization, stars with 
TD/D > 2 will then always have (independent on the value on 
Xjd) TD/D > 1 and can be regarded as thick disk stars. Stars 
with TD/D < 0.6 will always have (independent on the value 
on Xju) TD/D < 1 and can be regarded as thin disk stars. Stars 
having TD/D ratios in between these values will alter their 
classification as thin or thick disk stars as the normalization 
increase or decrease (in the range 6-14 %). We have therefore 
referred to these stars as "transition objects" throughout the pa- 
per. 
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